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Abstract 

giant quadrupole resonances are studied with the help of the Wigner Function Moments 
method generalized to take into account spin degrees of freedom. Equations of motion for 
collective variables are derived on the basis of Time Dependent Hartree-Fock equations in 
the model of harmonic oscillator with spin orbital mean field potential plus quadrupole- 
q ■ quadrupole residual interaction. Introducing spin allows one to consider new types of 

nuclear collective motion where the nucleons with spin 'up' oscillate against nucleons 
with spin 'down'. 
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The nuclear scissors mode was predicted [1, 2, 3, 4] as the collective motion of two types of 
nucleons - the protons undergo counter-rotational vibrations with respect to the neutrons. 
However, its collectivity turned out so small that it was even discussed whether [5] this mode 
is really collective. Pure phenomenological models (such as, for example, two rotors model [6]) 
could not help clarifying this question. The results of RPA calculations [5] were in qualitative 
agreement with the experiment, an indication in favour of non collectivity. The final conclusion 
was that the scissors mode is [7] "weakly collective, but strong on the single-particle scale". 
As a result [7]: "The weekly collective scissors mode excitation has become an ideal test of 
models - especially microscopic models - of nuclear vibrations. Most models are usually cali- 
brated to reproduce properties of strongly collective excitations (e.g. of J 71 " = 2 + or 3 _ states, 
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giant resonances, ...). Weekly-collective phenomena, however, force the models to make gen- 
uine predictions and the fact that the transitions in question are strong on the single-particle 
scale makes it impossible to dismiss failures as a mere detail, especially in the light of the 
overwhelming experimental evidence for them in many nuclei [8, 9]." 

The Wigner Function Moments (WFM) or phase space moments method turns out to be 
very useful in this situation. On the one hand it is a pure microscopic method, because it 
is based on the Time Dependent Hartree-Fock (TDHF) equation. On the other hand the 
method works with average values (moments) of operators which have a direct relation to the 
considered phenomenon. That makes it an ideal instrument to describe the basic characteristics 
(energies and excitation probabilities) of collective excitations such as, in particular, the scissors 
mode. Our investigations have shown that already the minimal set of collective variables, i.e. 
phase space moments up to quadratic order, is sufficient to reproduce the most important 
properties of the scissors mode: its inevitable coexistence with the Iso Vector Giant Quadrupole 
Resonance (IVGQR) implying a deformation of the Fermi surface. Also there is no doubt about 
its collectiveness, because the principal variables generating this mode (angular momenta of 
protons and neutrons) are collective being the respective operators averaged over all protons 
and neutrons. 

Further developments of the Wigner Function Moments method, namely, the switch from 
TDHF to TDHFB equations, i.e. taking into account pair correlations, allowed us to improve 
considerably the quantitative description of the scissors mode [10, 11]: for rare earth nuclei the 
energies are reproduced with ~ 10% accuracy and B(M1) factors were reduced about two times. 
However, they remain about two times too high with respect to experiment. The reason of the 
last discrepancy is probably hidden in the spin degrees of freedom, which were so far ignored 
by WFM method. One can not exclude, that due to spin dependent interactions some part of 
the force of Ml transitions is shifted to the energy region of 5-10 MeV, where 1 + resonance of 
the spin nature is observed. The generalization of the WFM method to take into account spin 
degrees of freedom is the goal of this paper. In a first step, we include in the consideration 
only the spin orbital interaction, as the most important one among all possible spin dependent 
interactions because it enters into the mean field. This allows us to understand the structure 
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of necessary modifications of the method avoiding too cumbersome calculations. In this way it 
becomes clear already on the level of formulation of the equations of motion for new collective 
variables, that we are faced with a new type of collective motion, namely the spin scissors mode. 
It turns out that the experimentally observed group of peaks in the energy interval 2-4 MeV 
corresponds to two different types of motion: "standard" scissors and this new kind of mode, 
i.e. spin scissors. 

The paper is organized as follows. In Sec. 2 the TDHF equations for 4-component density 
matrix are formulated and their Wigner transform is found. In Sec. 3 the spin structure of the 
density matrix and the model Hamiltonian are studied. In Sec. 4 the collective variables are 
defined and the respective dynamical equations are derived. In Sec. 5 the procedure of calcula- 
tion of excitation probabilities is recalled and the results of calculations of energies, B(M1) and 
B(E2) factors are discussed. Last remarks and the outlook are given in the conclusion section. 

2 TDHF equation with spin 

In this section we will consider the TDHF equation in coordinate space keeping all the spin 
indices. The TDHF equation in operator formulation reads [13] 

ihp=[h,p] (1) 

or in matrix form 

iTi < r, s|p|r", s" >= 

I d 3 r'(< r, s\h\r', s' >< r', s'|p|r", s" > - < r, s|p|r', s' >< r', s'\h\r", s" >). (2) 

s' J 

We do not specify the isospin indices in order to make the formulae more transparent. They 
will be re-introduced at the end. The set of TDHF equations (2) with specified spin indices 
reads 

ih < r\'p\r" > n = J d 3 r'(< r\h\r' > n < r'|p|r" > n - < r|p|r' > n < r'\h\r" > n 

+ < r|%' > n < r'\ P \r" - < r\p\r' > n < r'\h\r" 
ih < r\'p\r" > n = J d 3 r'(< r\h\r' > n < r'|p|r" > n - < r\p\r' > n < r'\h\r" > n 
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with the conventional notation 

| for s — ^ and | for s = — ^. 

These equations will be solved by the method of phase space (or Wigner function) moments. 
To this end we will rewrite the expression (3) with the help of Wigner transformation [13]. 

2.1 Wigner transformation 

The relevant mathematical details can be found in Appendix. To make the formulae more 
transparent we will not write out the coordinate dependence (r, p) of the functions. The 
Wigner transform of (3) can be written as 

ift/tt = ih {h^,f^} + h n f^-f n h^ + -{h^,f^}--{f\h^} 

-y{{^,/ 4t }} + y{{/ n ,^ t }} + -, 

ihp± = rHh n -h^)+^{(h^+h%f^}-^{{(h^-h^)j^}} 

-h n (f n ~ f U ) + y {h n , (f n + f u )} + y{{^ n , (/ n - f u )}} + .... (4) 

where the functions h, f are the Wigner transforms of h, p respectively, {/, g} is the Poisson 
bracket of the functions /(r, p) and g(r,p) and {{/, g}} is their double Poisson bracket; the 
dots stand for terms proportional to higher powers of h. Two more equations are obtained by 
the obvious change of arrows t^4- 

It is useful to rewrite the above equations in terms of functions / + = f^+f^, f~ = f^—f^- 
By analogy with isoscalar f n + f p and isovector f n — f p functions one can name the functions 
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/ + and / as spin-scalar and spin-vector ones respectively. We have: 

thj+ = j{h + ,f + }+^{h-,f-}+ih{h^j d }+ih{h^,r} + ..., 
in/- = j{h\f-} + l ^{h-,n-2h^r+2h^f d 

Mr = -h^r + h-r+ l ^{h^j + } + ^{h + j u } 

+y{{^r}}-y{{^,r}} + -, 
ihf d = h^r-h-f d +^{h^,n+^{h\f d } 

-y{{^r}} + y{{^,/ d }} + -, (5) 
where f u = f n , f d = f^, h± = h n ± h u . 



3 Model Hamiltonian 

The microscopic Hamiltonian of the model, harmonic oscillator with spin orbital interaction 
plus separable quadrupole-quadrupole residual interaction is given by 

A -2 -1 2 z N 

H = E(f^ + 2^-^)1 A] + r E (-irEEU'iW';) 

i=l /i=— 2 i j 

Y 2 Z N 

+ 2 K E (-!)"{£ ^( r *)^( r i) + E^r*)^)}, (6) 

^=-2 i^j iytj 

where N and Z are the numbers of neutrons and protons, respectively. The quadrupole operator 
q 2fl can be written as the tensor product: <Z 2 //( r ) = \Jl6n/5 r 2 F 2 /t(^, 0) = \/6{r <S> ?"}2,u, where 



r_i,r ,ri are cyclic coordinates [14] and C\^ Xv is the Clebsch-Gordan coefficient. The mean 
field potential for protons (or neutrons) is 

V T = \muj 2 r 2 + J2(-irZ^{r ® r} 2M - ^(r)iS. (7) 
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Here 

= /4P,r}r^r + (r,p,t) (8) 

with fd{p,r} = {27iK)- 3 jd 3 p jd 3 r. 
3.1 Spin structure 

Let us clarify the spin structure of the density matrix and the Hamiltonian. 

3.1.1 Density matrix 

As an example we consider the spherical case: 

^n«j m (r,s) =< nljm\r,s >= ft n ;(r) £C/™ 1ct Y;a(0, 0)xi CT (s), 



A,<T 

*>2" """ 2 L 



< Zjm (r, S ) =< r,a|nZjm >= ^(r) £ C^^ A (MM>)- (9) 

A,<7 



(10) 



Spin functions are defined as 

Xii ( s) = s XtW , xl . iW = (£) = x , (s) , X L( S ) = (f M , ). 

The density matrix reads [13] 

< r, s\p\r', s' > — £ v^ij < r, s\nljm >< nljm\r', s' > 

nljm 

= E vl l3 K nl {r)K nl {r>) £ C^^C*^ 

nljm ci,(T2 



where i;^ - are the occupation numbers of level nlj. According to the definition (10) the spin 
structure of the density matrix is 

«"^'>-«-">(K; ££)• (12) 

where 

p(rs, rV) = J2 v lijKijrrSj, s)(j)* nljm (r' , s') (13) 

nljm 
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and 

(f>nijm(r,s) = K n i(r)C 3 ™_ s j_Y lm _ s (6,4>). 



3.1.2 1-s Hamiltonian 

Written in cyclic coordinates, the spin orbit part of the Hamiltonian reads 
h ls = -r)(r) E (-) V-„ = - V (t) i 

where [14] 



] h t h 

m=-i V _< i75 ~ *02 



ii = ft(r Vi - riV ) = "^[ xV ^ _ + _ ^ V ?/)] = ~~^$ x + ^w)' 

/ = ft(r_iVi - riV_i) = -iU(xV y - yV x ) = l z , 

ft 1 

L 1 = ^(r_xV - r V_i) = ^[iV, - 2V3 - i(yV z - zV w )] = -^(4 - iZ„), 

Sl ~'7^{o o)' 5o "2 v o-iJ' ^"Tlli oj' 

The matrix elements of hi s in the configuration space are: 

< u\hi s \u' >= E J d 3 r < i/jr, s > hi s < r, s'|i/ > 

s,s' 

= -E /^WE^^l(M)xJ» 

77(r)E(-)^( r )^-M^'W E q^vWMk^M, 

A 4 A',<t' ' 2 

where v = n,l,j,m. Taking into account, that 

ft 

X{(s)5 1 xt(s') = 0, Xt( s )^-iXt( s/ ) = °> x\(s)S Xt{s') = ^^'i' 



xlW&X^) = 0, xi(s)£-iX+(s') = 0, xl(a)5 x;(s') = -^_i^_i, 



ft 



we find: 



< v\h ls \v' >= - Wd 3 r £ n^Y^e^Hr 

s,s' J A, A' 



/ / ft- 

iiv, 22 '22 Zi 2 2 



ft 



ft 



ft ; 
—i 

2 ■ ' ■ 2 2 2 

'ft: 

_,?7(rj 

2" 



ui, 2 2 ' ' 2 2 X 2 



TZ n >v{r)Y v 



ft 



E^(r)WM)c£f is , 



^o(r)5 s i<J 8 ,i + -^f_i(r)5 s i<J 8 ,_i 



ft - ft- 
-^=Zi(r)<J a _i£ 8 ,i - -Z (r)5 s _i^_i 



^n'l'j'm'ir, S f ) 

=< u\h n \u' > + < z/|ft n |z/' > + < u\h lt \u' > + < z/|ft w |z/ / >, (17) 



where 

ft - ft - 

ft n (r, s, s') = --?7(r)Zo(r)5 s i^i, ft n (r, s, s') = -^f|(r)Li(r)^i^i, 

Matrix elements of ft^ s in coordinate space are derived using the completeness of the set of 
functions (fi n ijm( r , s )- For example 



< ri, si|ft n |r 2 , s 2 >= E < r i> s iW >< v\h n \v' >< i/|r 2 , s 2 > 



E E 0n/jm(ri,Sl)E / rf3r C/jm( r > S )^W^o(r)5 st ^'t0n'«'/mK r > S OC'/'/m'( r 2,S2) 
nijm n'l'j'm' Q 



= -- Ey rf3 r<5(r - ri)5 SlS ?7(r)/o(r)5 st 5 s ' t 5(r - r 2 )5 S2S > = --r](ri)/ (ri)5(ri - r 2 )6\, lt <5 S2 t 



= ^/ (r 2 )h(r 2 )5(r 2 - ri)]«$ ait £ aat = /ij (ri,r 2 )<W<W, ( 18 ) 
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where 



Mri,r 2 ) = -^(rOZoCrO^n - r 2 ) = |f (r 2 ) [ V (v 2 )S(v 2 - n)]. (19) 



Analogously we find 

< ri,Si|/^|r 2 ,S2 >= -/ii (ri,r 2 )<5 ai 4.<J a2 4., 
< r 1 ,s 1 \hl^\r 2 ,s 2 >= hi_ 1 (r 1 ,r 2 )5 Slt 5 S2h 

< r 1 ,s 1 |^ t |r 2 ,s 2 >= h h {r 1 ,r 2 )5 Sll 5 S2t (20) 

with 

^ii(ri,r 2 ) = -^77(ri)fi(ri)<J(ri -r 2 ), h^r^r^ = --^77(ri)Z_i(ri)<S(ri - r 2 ). (21) 
As a result 

< r 1: s^htslr^ s 2 >= h lo (r 1: r 2 ) [£ Slt 5 S2t - 5 Sli 5 S2i ] + fi^fa, r 2 )5 Slt 5 S2 ± + h h (r u r 2 )5 Sli 5 S2t .(22) 

According to (52)-(54) its Wigner transformation reads: 

h ts {r, p) = -|?7(r){Zo(r, p)[5 Slt 5 S2t - 5 Sli S S2i ] + >/2Z-i(r, p)<S ait <W - V^h{r, p)5 sli 5 S2t }. (23) 

Taking into account equations (19) and (21) one can show that second and third equations 
of (3) are complex conjugated if the following relation is fulfilled: 

< r\p\r" > n =< r"\p\r > lt *, 

i.e. if the density matrix is Hermitian as it must be. 

4 Equations of motion 

Integrating the set of equations (5) over the phase space with the weights 

W = {r <g> p} Xfl , {r <g> r} AM , {p <g> p} v and 1 (24) 

one gets the dynamical equations for the following collective variables: 

L%(t) = Jd(p,v){r®p} x ,f ± (r,p,t), 
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Ll»(t) = Jd(p,T){r®p} x ,r 

R%(t) = Jd(p,v){r®r} x »f 
Rx»(t) = /d(p,r){r®r}^/» 

= 

Piit) - 



d(p,r){p®p} x J ± 
d(p,r){p®p} x J u 
d(p,r){p®p}xvf 
F ± (t) = fd(p,r)f- 
F u (t) = Jd(p,v)f 
F d (t) = Jd(p,r)f 



(r,P,t 
(r,P,t 
(r,P,t 
(r,p,t 
(r,P,t 

(r,P,t 
(r,P,t 

(r,P,i 
(r,P,t 
(r,P,t 
(r,P,t 



(25) 



We already named functions / + (r,p,t) and /~(r, p, t) as spin-scalar and spin- vector ones. So 
it will be natural to name the respective collective variables X£ (t) and (t) as spin-scalar 
and spin- vector variables. The required expressions for hr* 1 , and are 

h + = mu 2 r 2 + 2 £(-l)"Zg-{r ® r h-», h~ = -hrjl , h n = ~^=vU, h lt = ^=vh- (26) 

The integration yields: 

Lt = -P x + ,-mu 2 Rl + 2V5j2^2^{^ l }{Z^Rl}x, 

171 3=0 

-ih%L^ + ^-^(A + ^+l)^ + ^(A + / i)(A-/x + l)Lt_ 1 ], 



m 



P A - - mcu 2 ^ + 2^5 £ V 2 ^' + HSJH^" ® ^7 W 



j=0 



+irjV2 jd{p,v){rp) x ,[hr + Uf d ] 



TU 

^A/i+1 



—2 



V2 



jd(p,r)(rp) x , +1 [Uf - -V2l p 
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^t-i = -i 3 AVi- mw X,-i + 2v^Ey27TT{^}{z 2 + ®^} v _ 1 

_^| v /( A + ^)(A- / , + l)Lj 1 + -rjS^oF- - v^F**], 
F" = 2r ) [L d l _ l + L« n ], 

F d = -77^ + ^0], 
^ " m A ^ 

-ift^i^ + ^-^(A + ^+l)^ + v /(A + /i)(A- /U + l) J R^ 1 ], 
K = ^ L a, - + ir)V2 J d(p, rXr 2 )^/" + 

= ^Vi-^/ rf (p. r )(^)A, + i[Lir-v / 2/ r] 
^t-i = ^t-i-^/ rf (P' r )(^ 2 )A,-i[^r + v / 2^l 

_^| v /( A + ^)(A- / , + l)it:+ 1 , 

P+ = -2m W 2 L^ + 4v / 5E v /27 + T{^}{^ + ®L+} v 

i=o 

-ift^p- + y/iX-MX + n + l)^ + v /(A + /U )(A- /U + l)Pi_ 1 ], 

P A - = -2m^ + 4v / 5EV / 27 + T{2l{}{^2 + ®^7}A, 

i=o 

-^|^a + m + ivy/2 J dip, r)(p 2 ) x ,[hr + Uf d ], 

= -2m W 2 L^ +1 + 4v / 5Ey27TT{^}{Z 2 + ®LJ} v+1 

i=o 

^| v /(A-/i)(A + /i+l)P+ - ^ / d(p, r)(p 2 ) AM+1 [/_ 1 /- - v^of], 
P^_i = -2m W 2 Lt_ 1 + 4v^Ev / 2jTI{^}{Z 2 + ®^} v _ 1 
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-^(A + / .)(A-// + 1)P+ -i±J d(p, v){p 2 ) x ^[hr + V2l f% (27) 

where {l x {} is the Wigner 6j-symbol. For the sake of simplicity the time dependence of tensors 
is not written out. It is easy to see that this set of equations for moments of zero and second 
orders is not closed due to the integral terms of the type / d(p, r)(p 2 ) x ^i[lif~ + \^2l f d ] (see last 
equation). These terms generate moments of fourth order. They will be neglected according to 
the rules of the WFM method [15] . 

We are interested in the scissors mode with quantum number K n = 1 + . Therefore we only 
need the part of equations (27) with /i — 1. These equations are nonlinear and will be solved 
in small amplitude approximation. 

4.1 Linearized equations (fi = 1) 

Writing all variables as a sum of their equilibrium value plus a small deviation 

Rx»(t)=RZ + n x »(t), P Xlx (t) = P e x l + V x »(t), L x ,(t) = Lll + C x ,(t) 

and neglecting quadratic deviations, one obtains the desired equations. Naturally one needs 
to know the equilibrium values of all variables. Obvious equilibrium conditions for an axially 
symmetric nucleus are: 

Rt ±1 (eq) = Rt± 2 (eq) = 0, Rt (eq) ^ 0. (28) 

It is obvious that all ground state properties of the system of spin up nucleons are identical to 
that of the system of nucleons with spin down. Therefore 

RxM) = P A " (eg) = L v (eg) = 0. (29) 

We also will suppose 

L^eq) = L^(eq) = L^eq) = and R^(eq) = R%{eq) = 0. (30) 
With the help of the above equilibrium relations one arrives at the following set of equations: 
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-^|[£ 2 - 1 + 2£« 2 + 2^ ], 

-^P 2 - - mu 2 n 21 - ^Z+(eq)TZ 21 - ih^C^, 

^V u 22 -mu 2 n u 22 - sj\z 2 \(eq)K u 22 - ih^£+, 

±V d -mu 2 TZ d 2() - ^Z+(eq)TZ d + ^=Z^eq)TZ d m - ih\ 

^-[RUe q )Z+ - Z^eq)H\ x \ - ih^ + V2£?„], 

-\Hz 2 + (e q )n 21 - ih^Cf, - ^- V F d , 

-T^n + V^o], 

-c 21 - ih^n^, 
— ry- -jh-TZ + 

*~22 <2 21 ' 



2 

m 



^20 2V2 

- ih 7 L [V - + 2V2 2 + 2^V d }, 

-2mu 2 C 21 + VGZ^eq)^ - ^Z+(eq)£ 21 - 

-2ma; 2 ^ 2 + 2^Z^eq)C u 22 - i7i|7>+, 

-2m^ - 2y|z 2 + (eg)^ + 2^|z 2 + (eg)^ - i^^f* 



^P d - mc^Xo + ^^2 + o(eg)^ , 
-2mw 2 £g + 2 J^Z+(eq)£ d 20 . 
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Let us recall that all variables and equilibrium quantities R£o(eq) and Z 20 (eq) in (31) have 
isospin indices r = n, p. All the difference between neutron and proton systems lies in the 
mean field quantity Z 2 Q~(eq), which is different for neutrons and protons (see (8)). 



4.2 Isovector, isoscalar 

It is convenient to rewrite the equations (31) in terms of isovector and isoscalar variables 

p pn pP p pn pP f rn rp 

-«A M - -Ka m _ A/i' - -*A/i ~ ^A/x> ^A m - -^A M _ ^A/x" 

It is natural to define also isovector and isoscalar strength constants Xi = \{x~ X) an d Xo — 
\{x + X) connected by the relation xi — a Xo [16]. Then the equations for the neutron and 
proton systems are transformed into isovector and isoscalar ones. The equations for the isovector 
system are given by 



-21 



v+ - 

• 21 



m 



2 1 
mu 2 - -j=axoRo (eq) + -^=(1 + a)xoRto(eq) 



/v 21 



r — 

•-U 

r — 

■-d 

£■20 = 

Ai = 

Ki = 



F 

^ 4 

n, 
n 



21 



- ih 7 l[jC 21 + 2jC, u 22 + 2 ] j-jC d 20 }, 

1 1 T) 

—V 21 - [mu 2 + — =x R2o(eq)]Tl 21 - ik-C^i, 
m a/6 2 



' 22 

m 



2 - 71 - 

[mu 2 - \ -XoR2o(eq)]Tl 22 - ih-C^, 



±-V d -\muj 2 + ^XoR2 (eq)]n d 20 + ^=xMeq)K d m - ih^^, 

-y|(l - a)xoR2o(eq)n+ 1 - ih^[C^ + V2C d w ], 

-^lxoR2o(eq)'Jl 21 - h^[iC+ + hF*\, 
V 



C 10 = -h-^liCi, + hF% 



21 



'2a/2 
-^n + V^Cfo], 

l£+ _ + 27^ 2 + 2^|^ ], 

Z_ 2 1 " i '^-21) 

m 2 



14 



— u 2 - T) - , 

^"22 = ~^22 ~ ^2^ 21 ' 

^2i = -^mc^ + ^Xo^oMK^ + v^^ 

^_ _ n 

P 21 = -\2mu 2 + \l -XoR2o(eq)}£ 2 ~ 1 + \^XoR2o(eq)C u -ih-V^, 

P 22 = -[2muj 2 -2^ Xo R 20 (eq)}Q 2 -ih^V 2 + 1 , 

P 20 = -[2mu 2 + 2^XoR2o(eq)]Ci + ^ X oR2o(eq)^ - l h^P 2 \, 

' d \ 2 

A)o = -Poo-^^oo + ^Xoi?2o(eg)^ , 

■-d 2 - A 

00 ~~ m 00 ' 

lit 

P d 00 = -2mu 2 C d 00 + ^= Xo R 20 (eq)C d 20 . (32) 
The isoscalar set of equations is easily obtained from (32) by taking a — 1. 



4.3 Angular momentum conservation 

The set of equations (32) contains three integrals of motion (see Appendix B). The first one is 
(in the isoscalar case) the conservation of total angular momentum < J ± >—< l ± > + < S± >. 
By definition 

< h >= Tr(hp) =J2 j d?>r j d * r ' < r \hW >< r, s|p|r, s > 
= J2 J ^ r j rfVfi(r)5(r - r') < r', s\p\r, s >= J d 3 rh(r)[< r\p\r > n + < r|p|r > u ] 

= J d(pr)h(r, p)/ + (r, p) = -iV2 J rf(pr){r ® p} n / + (r, p) = -iV2L+. (33) 

The average value of the spin operator Si reads: 

< 5i >= TV(Sip) = J2jd 3 r< soils' >< r, s'|p|r, s > 

s,s' 

= E< > | rf( P r)r' s (r,p) = — ^ st ^ = = _JL F <* (34) 
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As a result < J\ >= —-^(2iLi 1 -\-hF d ). It is easy to see that such combination of the respective 
equations in (32) is equal to zero in the isoscalar case (a = 1), i.e. the total angular momentum 
is conserved. 



5 Energies and excitation probabilities 

Imposing the time evolution via e lVit for all variables one transforms (32) into a set of algebraic 
equations. Eigenfrequencies are found as the zeros of its secular equation. Excitation proba- 
bilities are calculated with the help of the theory of linear response of the system to a weak 
external field 

6{t) = 6 exp(-iQt) + 6 f exp(iQt). (35) 

The detailed explanation can be found in [16]. Here we only will recall the main points. The 
matrix elements of the operator O obey the relationship 



< ij a \O\ij > | 2 = h lim (ft - Q a )< ip>\0\ip' > exp(-ifti), 



(36) 



where ipo and ip a are the stationary wave functions of the unperturbed ground and excited 
states; ip' is the wave function of the perturbed ground state, Q a = (E a — Eq)/% are the normal 
frequencies, the bar means averaging over a time interval much larger than 1/fl. 

To calculate the magnetic transition probability, it is necessary to excite the system by 
the following external field: 



CV = A*jv (g 8 &/h - Wiy^i* x V]) V(r A F v ), fi N = 



(37) 



Here gi — 1, g s — 5.5856 for protons and g t — 0, g s — —3.8263 for neutrons. We are interested 
in the dipole operator (A = 1, // = 1). In cyclic coordinates it looks like 




gA/h-giV^^Cll^V, 



(38) 



Its Wigner transformation is 



(Ou)w — fJ'N 




g s S! - iV2~ gi CH M r uPa 



/ft, 
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For the matrix element we have 



< ^|6ii|V' >= fiN\ — [<7«<^l W > -iV2L p +] /h 

47T L J 



.?-(g»F d » + g»F d P)-iLl+ 



/% 



UN 




-\M + 9l)F d + (rf - gl)F d ] - - £+) 



^G? S P - ^ d + ^1 + ^ + - 1]^ 



(39) 



Deriving (39) we have used equilibrium relations (30) and 2iL\ x = —TiF d , which follows from 
the angular momentum conservation. 

Due to the external field one of the dynamical equations of (32) becomes inhomogeneous: 



TZ 21 - -£+ + %h\\RT 2X + 2K U 22 + 2d^n d 20 ] = t^p-R 2Q (eq) exp 



N 



JQt 



m 



(40) 



Solving the inhomogeneous set of equations (32) in the isovector and isoscalar cases one can 
find the required values of Cfi , Cti and F d and calculate B(M1) factors for all (isovector and 
isoscalar) excitations. 

To calculate the electric transition probability, it is necessary to excite the system with 
the external field operator 

O v = er 2 r v = P{r ® r} 2 ^, (41) 

where (3 = ey^. Its Wigner transform is identical to (41): (0 2 n>)w = f3{r®r} 2 ^. We consider 
the case // = 1. The matrix element is given by 



< >= mil = \mii - Rii) = \mti - ^ii)- 



(42) 



The external field yields two inhomogeneous equations of (32): 



^21 ™^21 + 

m 



moo 



-^axoRooieq) + -Wl + a)xoRl {eq) 



21 



Mhl\C 21 + 2C u 22 + 2 x r-C. 



20J 



'11 



^3 (x i2w(e9) " R °° {eq) ) 6XP ^' 
\{l - a) X oR2o(eq)nl 1 + ih V -[C u + V2C d l0 ] = -^R 20 (eq) exp 4 ™ 



(43) 
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Solving the inhomogeneous set of equations (32) in the isovector and isoscalar cases one can find 
the required values of TZ^i and H21 and calculate B(E2) factors for all (isovector and isoscalar) 
excitations. 



5.1 Isovector excitations 

The energies and excitation probabilities obtained by the solution of the isovector set of equa- 
tions (32) are given in Table 1. These results will be discussed by comparing them with some 
simplified cases, that will allows one to understand the physical nature of the considered exci- 
tations. 

Tables 1, 2. Energies and excitation probabilities of the isovector system, calculated for 
164 Er with two values of the spin orbit interaction constant rj. The quantum numbers (including 
indices +, — , u, d) of variables responsible for the generation of the present level are shown in 
a first column of Table 2. 



Table 1. rj = 0.361 MeV 



E iv , MeV 


B(M1), nl 


B(E2), B w 


(A,/i) s 


E iv , MeV 


B(M1), nl 


B(E2), B w 


1.61 


3.54 


0.12 


(1,1)- 


1.70 








2.18 


5.33 


1.02 


(1,1) + 


2.07 


9.03 


1.17 


12.80 


0.01 


0.04 


(0,0) d 


12.81 








14.50 


0.01 


0.03 


(2,1)- 


14.50 








16.18 


0.02 


0.18 


(2,2)" 


16.20 








16.20 








(2,0) d 


16.20 








20.59 


2.78 


35.45 


(2,1) + 


20.55 


2.81 


35.70 


0.26 i 


-5.43 i 


0.0 i 


(l,0) d 










Table 2. rj = 



5.1.1 The limit rj = 



In the case of vanishing spin orbit potential, 77 = 0, the set of equations (32) splits into four 
independent subsets. 
The first one is 



-21 



1 21 

m 



1 



mu 7 =«Xo-Roo(eg) + —7= (I + a)xo-R2o(e<?) 



^3 



/V 2H 
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Ai = -y |(! - a)xo-R2o(eg)^2i, 

= -[2mw 2 + ^|xo^2o(eg)]£++v / 6xoi?2o(eg)^ 1 . (44) 

This set of equations for spin-scalar variables coincides with the set of equations (29) of [16] 
and describes the joint dynamics of the "standard" nuclear scissors mode and the /i — 1 branch 
of IVGQR. Remembering that Xo — 6k , R 2fM = Q2^/V^, Q20 = Qoof^ an d taking the self 
consistent value k = — f-^ [17] we find the energies of these modes 

Q% = u 2 {2 - a)(l + 5/3) ± y/u*(2 - a) 2 (l + 5/3) 2 - 4w 4 (l - a)5 2 . (45) 

Usually [16] we choose a = —2. Then 

nl=4tf 2 |l + ;±t/(l + J)'-^|. (46) 
The second subset reads 




r — 

J - / 21 ~ 


!^ r 2 1 

— "r 2] - raw H — =■ 


Ai = 


-]J^XoR2o(eq)n 21 , 






C = 


lit 




m 


K = 


-|2mw 2 + y jjxo#20 



'11- 



(47) 



In the absence of spin dependent forces one could naively expect that this set of equations 
should be identical with (44). However we observe an essential difference between the two, the 
difference being determined by the variation of the mean field (the variable Z 2 i)- This fact can 
easily be understood by comparing the first and second equations of (27). The terms responsible 
for the contribution of the mean field are {Z 2 <g> R 2 and {Z 2 <8> R 2 }\^- Their variations are 
{SZ^ <g> Rt(eq) + Z^(eq) <g> 5R%} Xtl and {5Z£ <g> R 2 (eq) + Z^{eq) <g> 5i?2 }a m respectively. We 
assumed that R 2fl {eq) = 0, so the contribution of the mean field variation 5Z 2/X in the second 
equation disappears and we arrive to the announced difference. In other words, when the 
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particles with different directions of spins move in phase, the mean field variations generated 
by their motion add and, when they move out of phase, the respective mean field variations 
annihilate one another. 

The set of equations for spin-vector variables (47) has the following eigenfrequencies: 



The numerical estimates for l64 Er are shown in Table 2. The low lying level has E^low) = 
%Vt_ = 1.70 MeV and the high lying level has E^ihigh) = hQ + = 14.50 MeV. So, in addition to 
the well known ("standard") scissors mode with energy E^ = 2.07 MeV, describing the relative 
rotational vibrations of all (spin up + spin down) neutrons with respect of all protons, we get 
one more scissors mode, which is generated by the new type of the nuclear collective motion - 
the relative rotational vibrations of spin up nucleons with respect of spin down nucleons. The 
inclusion of the spin orbital interaction only slightly shifts these levels apart (compare Tables 
1 and 2): 1.70 — > 1.61 and 2.07 — > 2.18. However, the main role of the spin orbital interaction 
consists in the excitation of the new scissors mode. This mode is not excited directly by the 
electromagnetic field (see formula (40)). It can be excited only indirectly, via the coupling with 
the "standard" scissors mode by means of the spin orbital interaction. 

The new scissors mode is generated by the variable responsible for the relative rotational 
motion. It is easy to see by removing this variable from the set of equations (47). In this case 



its characteristic equation gives only one level with energy E = 2huyl + 5/3, which practically 
coincides with E^high) and lies in between the isoscalar and isovector giant quadrupole res- 
onances. From the physical interpretation of variables ^-2i) ^2i> ATi > responsible for this mode, 
it follows that it describes rather complicate motion: the proton system oscillates out of phase 
with the neutron system, whereas inside of each system spin up nucleons oscillate out of phase 
with spin down nucleons. In the absence of the proper residual interaction it has the standard 
shell model value of the energy E = 2hu(5). So, the new scissors mode is accompanied by the 
high lying excitation, which can be called the spin-vector giant quadrupole resonance. As one 
can see, the situation is quite similar to that of the "standard" scissors mode which exists only 
together [16] with isovector giant quadrupole resonance. 




(48) 
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Let us point to the interesting feature of the new scissors mode. As one knows, the variable 
C u is spin- vector (see the text after eq. (25)), i.e. it describes the relative rotational oscillations 
of spin up nucleons with respect of spin down nucleons (spin scissors). This means that we have 
proton spin scissors and neutron spin scissors. At the same time is the isovector variable, 
i.e. it describes the relative rotational oscillations of protons with respect of neutrons. So, 
the resulting motion described by the variable L\ x looks rather complex - proton spin scissors 
vibrate in a rotational way (like scissors) with respect of neutron spin scissors, i.e. we have the 
scissors made of two other scissors! 

It is necessary also to stress the important role of the Fermi surface deformation (variable 
7^2i). Here the situation is exactly the same as in the case of the "standard" scissors mode. If 
one removes the variable P 21 from the set of equations (47), one gets a zero value for E^low) 



and a considerably changed expression for the high lying mode - E 2l {high) = hu)\j2(l + 5/3). 
The third subset reads 

K2 = ^P% 2 -[rnuj 2 -J^XoR2o(eq)}K 2 l 2 , 



2 

m 



7? — r u 

/v 22 — ''-'22) 



_ u 



V,, = -[2mco 2 -2J^ X oR2o(eq)]C u 22 . (49) 



Its solution gives E 22 = 2fnuy(l + -^S). 
The fourth subset reads 



Ao = ^^-[rnuJ 2 + ^XoR2o(eq)}ni + ^=xoR2 (eq)n d 00 , 



m 

d 2 



7? — r° 

20 — ™ 20) 



, d 

m 



V d 20 = -[2mu 2 + 2^xoR2o(eq)}Ci + ^=XoR2 (eq)C d 00 , 

d \ 2 

4o = -V d -mu 2 n d 00 + -j=xoR2o(eq)n d 20 , 

•- d 2 — j 

7? — — C d 

lit 

P d 00 = -2mu 2 C d + ^=XoR2 (eq)C d 20 . (50) 

Its solution gives E 20 = 2hui^J(l + |5) and E 00 = 2hui^J(l — p). 
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According to the physical interpretation of variables of the third and fourth subsets they 
describe the spin-flip modes: giant quadrupole and monopole spin-flip resonances with energies 
E 20 = E 22 = 16.20 MeV and E 00 = 12.81 MeV respectively (Table 2). Our Hamiltonian does 
not contain any forces, which can change the direction of the spin, so these modes have zero 
values of excitation probabilities and the standard shell model values of energies E = 2hu(S). 
The spin orbital interaction changes this result only negligibly. 

And finally, a few words about the solution with the imaginary eigenfrequency (last line of 
Table 1). From where does it appear? We carefully checked that it is not a mistake of calcu- 
lation or computation, because the energy weighted sum rule (see Appendix C) is fulfilled, the 
contribution of this unexpected mode being ~1.5%. What could be its physical interpretation? 
The time dependence of variables in this mode is e~ Qot (Q = 0.26 MeV), i.e. this excitation 
is unstable, i.e. decaying. This may be connected with the neglect by higher order moments. 
One can hope also that the inclusion of the spin-spin or spin-multipole residual interaction will 
produce the necessary restoring force and stabilize the mode. 

5.1.2 The approximation = = 

The set of equations for the variables X ± is obtained from (32) by setting variables X u ,X d 
to zero. Its solution gives energies and excitation probabilities given in the Table 3. The 
comparison with the exact solution shows that the approximation of two (spin up and spin 
down) "liquids" can be used with great care - it works very well to describe energies and B(E2) 
factors, but leads to erroneous results for B(M1) factors because the essential contribution into 
their values stems from the (now neglected) variable F d . 

Table 3. Energies and excitation probabilities of the isovector system, calculated for 164 Er 
with rj = 0.361 MeV in the approximation = = 0. 



E iv , MeV 


B(M1), pi 2 N 


B[E2), Byy 


1.62 


1.30 


0.16 


2.15 


7.74 


1.01 


14.50 





0.03 


20.59 


2.81 


35.66 
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5.2 Isoscalar case 

The energies and excitation probabilities obtained by the solution of the isoscalar set of equa- 
tions (32) are demonstrated in the Table 4. Analogously to the isovector case the results will 
be discussed in comparison with rj = limit, what helps to understand their physical nature. 

Tables 4, 5. Energies and excitation probabilities of the isoscalar system, calculated 
for 164 Er with two values of the spin orbital interaction constant rj. The quantum numbers 
(including indices +, — , u, d) of variables responsible for the generation of the present level are 
shown in a first column of Table 5. 



Table 4. 77 = 0.361 MeV 



E is , MeV 


B(M1), /4 


B{E2), B w 


(A,/z) s 


E is , MeV 


B(M1), fx% 


B(E2), B w 


1.73 


-0.07 


1.12 


(1,1)" 


1.70 








0.39 


0.24 


117.19 


(1,1) + 









12.83 





0.66 


(0,0) d 


12.81 








14.51 





0.12 


(2,1)- 


14.51 








16.20 








(2,2)" 


16.20 








16.22 





0.20 


(2,0) d 


16.20 








10.28 





66.50 


(2,1) + 


10.33 





67.47 


0.20 i 


-0.12 i 


30.30 i 


(i,o) d 










Table 5. rj = 



5.2.1 The limit of vanishing spin-orbit potential (r/ = 0) 

In the isoscalar case the set of equations (44) is transformed into 



^21 



Ml 



1 

m 
0, 



V+ - 

' 21 



mu 2 i=Xo^oo(eg) + ^=Xo#2o(e<?) 



/V 2H 



/V 21 — _' L '21, 



m 



^/|xo^2o(eg)]£ + 1 + v / 6xo^2o(eg)£ + 1 . 



(51) 



These equations coincide with the set of equations (24) of [16]. They describe the joint dynamics 
of the fi = 1 branch of ISGQR and the orbital angular momentum (its l\ projection). In the 
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absence of spin the orbital angular momentum is conserved < l\ >= —iy/^C^ = 0) as it 
should be. However, in the case with spin the situation is changed radically. Now the orbital 
angular momentum and spin can vibrate out of phase keeping their sum (the total angular 
momentum) unchanged. This new circumstance reveals itself in the appearance of the new 
low lying mode which describes the relative motion of the orbital angular momentum and 
spin of the nucleus (see Table 4). According to our calculations this mode (E=0.39 MeV) has 
essentially electric character with B(E2) = 117 W.u. demonstrating also a small admixture of 
the magnetic properties with B(M1) = 0.24fi 2 N . 

As it is seen, the set of equations (47) does not depend on a, so it is not changed in the 
isoscalar case. Naturally it has just the same eigenvalues: E(low) = 1.70 MeV, E(high) = 14.50 
MeV. The collective motion corresponding to these two modes is more simple, than in the 
isovector case. The variable £7i> responsible for the low lying excitation, describes the real spin 
scissors: all spin up nucleons (protons together with neutrons) oscillate rotationally out of phase 
with all spin down nucleons. This mode has B(M1)=0 because magnetic moments produced 
by two parts of spin scissors annihilate one another. The variables TZ^i , V21 , £21 > responsible 
for the high lying excitation, describe out of phase oscillations of all spin up nucleons (protons 
together with neutrons) with respect of all spin down nucleons. The inclusion of the spin 
orbital interaction influences practically only on the low-lying mode - its energy becomes a 
little bit higher (1.70 MeV ->■ 1.73 MeV) and it acquires the nonzero B(E2) value (B(E2)=1.12 
W.u.) demonstrating its electric character. The negligibly small negative B(M1) value and the 
imaginary solution (last line of the Table 4) appear, probably, due to the same reasons, as was 
already discussed at the end of section 5.1.1. 

5.2.2 The approximation f# = f n = 

Removing the variables X u ,X d from (32) with a = 1 we obtain the isoscalar set of equations 
for X ± variables. The results of calculations are shown in the Table 6. It is readily seen, that 
approximate results reproduce rather well the exact ones (Table 4) except for one case: instead 
of the excitation with the energy E=0.39 MeV, which is generated by the relative motion of 
the orbital angular momentum and the spin, we obtain an imaginary solution. This is not 
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surprising, since after removing the variable F d , one can not ensure the conservation of the 
total angular momentum (see section 4.3). 

Table 6. Energies and excitation probabilities of the isoscalar system, calculated for 164 Er 
with i] = 0.361 MeV in the approximation = = 0. 



E- iy , MeV 


B(M1), pt% 


B(E2), B w 


0.18i 


0.79 i 


-212.4 i 


1.72 


0.08 


0.99 


10.32 





67.34 


14.51 





0.12 



6 Conclusion 

In this work, the WFM method is applied for a first time to solve the TDHF equation including 
spin dynamics. The model Hamiltonian consists of a harmonic oscillator with spin orbit mean 
field potential plus quadrupole-quadrupole residual interaction. Spin dependent collective vari- 
ables are defined and a corresponding closed set of nonlinear dynamical equations is set up. 
The equations are solved in small amplitude approximation. Two isovector and two isoscalar 
low lying eigenfrequencies and five isovector and five isoscalar high lying eigenfrequencies have 
been found. Three low lying levels correspond to the excitation of a new kind of mode, un- 
known earlier. For example the isovector level with energy E=1.61 MeV describes the rotational 
oscillations of nucleons with the spin projection "up" with respect of nucleons with the spin 
projection "down", i.e. one can talk of a nuclear spin scissors mode. Thus the experimentally 
observed group of 1 + peaks in the interval 2-4 MeV, associated usually with the nuclear scis- 
sors mode, in reality consists of the excitations of the "spin" scissors mode together with the 
"standard" scissors mode. It is necessary to point out an interesting feature of this new mode. 
Being of isovector nature, it describes also the motion of protons with respect to neutrons and, 
therefore, the resulting motion can be characterized as the rotational oscillations of proton spin 
scissors with respect of neutron spin scissors. In the analogous isoscalar mode with the energy 
E=1.73 MeV the proton and neutron spin scissors move in phase. This mode has B(M1)=0 
because magnetic moments produced by two different parts of proton spin scissors annihilate 
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each other. A small B(E2) factor is explained by the coupling with the ISGQR. 

One more new low lying mode with energy E=0.39 MeV is generated by the out of phase 
oscillations of the orbital angular momentum and spin of the nucleus. The rather big B(E2) 
factor of this mode still waits for its explanation. 

There are ten high lying excitations, the two of them being really new: isovector and 
isoscalar " spin- vector" resonances with energies E=14.5 MeV. Further six high lying modes 
can be interpreted as spin-flip ones. In the absence of residual interactions all these eight 
modes have very small excitation probabilities and standard shell model energies E=2huj(8). 

In the light of the above results, the study of all discussed (low and high lying) excitations 
with proper residual interactions included will be the natural continuation of this work. It is 
known very well that pairing is very important for the correct description of the "standard" 
scissors. Therefore one has to take into account pair correlations too, i.e. to move from the 
TDHF equations to TDHFB equations. This also will be a task for the future. 
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7 Appendix A 

According to the definition [13] of Wigner transformation one finds 



Usually we will take r)(r) = rj(r). Then -£-r}{r) = y As a result (xS- — y-^-)r](r) = and 




(52) 



M r >P) = -^V(r)(xp y - yp x ). 



(53) 



In a similar way one finds 



M r >P) 



-r](r)[i(xp z -zp x )- (yp z 

h 

-r](r)[i(xp z - zp x ) + (yp z 



ZPy)}, 



fy-i(r,p) 



ZPy)}. 



(54) 
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The Wigner transformation of a product of two operators is given [13] by the following 
formula 



(hp) w = %,p)exp(yA)/(r,p) 

= h(r, p)/(r, p) + ^{h, /} - ^{{h, /}} + o(h 3 ), 



(55) 



where A=VrV P — V P Vr, \h,f} = hAf is the Poisson bracket of functions /i(r, p) and 
f(r, p), {{/i, /}} = h(r, p)(A) 2 /(r, p) is their double Poisson bracket. 



8 Appendix B 



Integrals of motion. 



v. 



d 
20 



-"21 



—2 



2(1 - a)xo . \^ cd . « ^ d 



ft ^ 
2 



hr]Rt (eq) V 2 



1 - 

77/ 



A" 2 + U^ 2 - ^XoRUeq) ) (^22 - C d ) 



« — ?nw 



2 2a X0 D+ , x . 

i?oo(eg) + 



(1 + a)xo 



V6 



Rt (eq) C a 



. hrjm 



3V2 



/2_, (i% + 
3 20 v / 2mXo^2o(e'?) 



where 



9 Appendix C 

The sum rule for magnetic transitions. 

Magnetic operator 

/ 3 

°h = \ ~r(9ih + QsS^HN/h (56) 

V 47T 

obeys the following sum rule 

- £ )(| < i/|On|0 > | 2 + | < i/|Oi-i|0 > | 2 ) = - < 0|[6 n , [//,6i_i]]|0 > . (57) 
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The oscillator part h osc of the Hamiltonian (6) commutes with S^, so we can use for the double 
commutator of with h osc the result found in [16]: 

[d 14> , [h osc , 6 14> ,}} = ^EEE {-ircH^C&^in ® r,} 2e {r, ® rj } 2 y N . (58) 

i j v,a,e 

Its ground state matrix element is [16] 

9 1 — a 

< 0| [On, [/w, Oi_i]]|0 >= ^X(^) 2 ^ = -^^Qoo^V/i^. (59) 

Now we have to calculate the double commutator of with the spin orbital part h\ s of the 
Hamiltonian (6). One has 

h ls = -v(r) E (-m,s_„ [/;, g = -V2C{; M U [s„ s^} = -V2Cl; M s u , 

M=-l 



1 1 

/i5 



z^=— i i 



[h ls ,St\ = -V(r) E = v^(r) E 

/Lt=— 1 A t=— 1 



fa., <V]] = -^J^(r) E(- 1 ) M <i^(^C'5,i^^^-m + ^C'5,i-m^^') 

= -3^^(r)(^-<? s )E{}i\^ (6°) 



\x=— 1 

= 3^y| ^(r)G7, - ^) EGiiX- 1 ® ( 61 ) 



[O,, 6^]] = -9^(r) 5>i - 5>i - ^(-lyiCii}^?,!^ ® # W ( 62 ) 

i — 1 
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We need the ground state matrix element of this double commutator for = 1, = — 1: 

<0|[6n,[fc|.,Oi-i]]|0> = -^tyi-9l)U9i-9l(-m]ll} 

Cf\ 11 <0|7 7 (r){r®^Lo|0>. (63) 

Averaging of spin variables implies the trace over spin indices, which is obviously equal to zero 
in our case. Therefore the spin orbital interaction does not contribute into the sum rule. 
The sum rule for electric transitions coincides with that of the paper [16]: 

J2(e„-e )(\ <u\6 21 \o> | 2 + | <v\6 2 -i\o> | 2 ) = -<o|[d 2 i,[if,6 2 _i]]|o> 

V 

= e 2 -^Qoo(l + 5/3). (64) 
m Air 
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